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Q ■ Abstract 

Exact solutions of the steady resistive three dimensional (3D) magnetohydrodynamics (MHD) 
i-£h ' equations in cylindrical coordinates for an incompressible plasma are presented. The solutions are 

translationally invariant along one direction and in general they describe a process of reconnective 
annihilation in a curved current layer with non vanishing magnetic field. In the derivation of 
the solutions the ideal case with vanishing resistivity and electric field is considered first and then 

c/3 ■ 

generalized to include the effects of finite electric field and resistivity. Particular attention is devoted 
£>~v to the analysis how the latter ones influence the presence of singularities in the solutions. In this 

respect comparisons with the purely two-dimensional case are made and the resulting important 

y—{ ' differences are highlighted. Finally, applications of the solutions for modeling an important class 

>'■ 

Q\ ■ of solar flares are discussed. 
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I. INTRODUCTION 



Magnetic reconnection is an effective mechanism for restructuring the magnetic field and 
converting magnetic energy into other forms of energy in plasmas with a high electric con- 
ductivity. This process plays a key role in large scale cosmic phenomena such as solar 
flares and geomagnetic substorms. Since the pioneering works of Dungey!, Parker!, Sweet! 
and Petschek!, where the basic physics of magnetic reconnection was clarified, many ex- 
act two-dimensional (2D) solutions describing some simplified versions of this process were 
found!!!!. 

Moreover, over the last few years there was a similar progress in the theory of three- 
dimensional (3D) reconnection. In particular, some simplified forms of this process at null 
points of 3D magnetic fields have been found as well §00. 

In this paper we present a generalization of exact solutions for a 2D curvilinear geometry!, 
the so-called two and a half dimensions (2|D) case, where a translational invariant compo- 
nent of velocity and magnetic field along the third direction exists. These solutions describe 
steady incompressible resistive magnetohydrodynamics (MHD) flows in a configuration with 
non-vanishing magnetic field. In some respects they resemble the 2D solutions describing a 
particular type of magnetic reconnection which is called reconnective annihilation!. So we 
will also use this term further on to denote the process described by our solutions. 

It should be noted also that when the present article was in preparation, a paper of 
Watson and CraigEl appeared, where similar solutions have been presented. Since we found 
these solutions independently, our considerations and interpretations differ in many respects. 

In sections [H] and |ITT| , respectively, the basic equations and the form of solutions are 



described. In section IV we derive the solutions in the limit of ideal MHD and discuss their 



properties, while in section [V| we consider how these solutions are modified in the case of 



non-vanishing resistivity. The conclusions are drawn in section VI 



II. BASIC EQUATIONS 

The set of MHD equations for an incompressible plasma with uniform density and resis- 
tivity consists of the equation of motion 

(v V)v = -Vp+ (V x B) x B, (1) 
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the Ohm's law 

E + vxB = i)VxB (2) 
and the divergence-free conditions for the velocity v and the magnetic field B 

V-v = 0, 

(3) 

V-B = 0. 

All equations here are written in a dimensionless form such that B and v are normalized 
to B e and v^ e respectively, which are characteristic values of the magnetic field and of the 
Alfven velocity. E represents the dimensionless electric field while r\ corresponds to the 
inverse magnetic Reynolds number. 

The current density j is determined separately by Ampere's law 

j = V x B. (4) 

Consider a cylindrical coordinate system (r, 9, z) where r and 9 are related to the Cartesian 
coordinates (x, y) in the following way 

x = rsin9, y = rcos9 — d, (5) 

where d > so that the pole is below the plane y = 0. Assume that in this coordinate 
system the functions B and v depend only on r and 9. Then they can be written as follows 

„ ^ ( ldA dA rr\ , \ fldi/j dip \ 

where A, H, ip and V are functions of r and 9 which are to be found. In particular, for 
H = V = this is the usual representation of two-dimensional magnetic and incompressible 
velocity fields in terms of a flux function A and a stream function ip, respectively. 

By using this representation we obtain from Eq. ([[]) that the functions A, H, ip and V 
must satisfy the following two equations 

[ij,V 2 ij] = [A,V 2 A], (7) 

[V, - [H, A] = Cr, (8) 

where C is an arbitrary constant and the Poisson brackets are used such that 

r, i = 9f_9g _ dg_d£ 
[1,91 8rd9 drd9' 
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Similarly, Eq. (0) gives the following two equations 

+ [A, V] =r]rV 2 H, 



(9) 



E^r + [ip, A] = -r]rV 2 A. 



(10) 



III. FORM OF THE SOLUTIONS 



For the system (m)— (10) we seek solutions of the form 



A(r,9)=A 1 (r)6 + A (r), 
ip(r,9)=ipx{r)9 + ip (r), 
H(r,9) = Hx(r)9 + H (r), 
V{r,6) = V 1 (r)9 + V (r). 

This form is a generalization of the ansatz used in Ref. ||| for a two-dimensional configuration. 
Substituting ([□]) into ©-(0) provides four equations each of which is a polynomial linear in 
9. Thus for each equation the part of the polynomial not depending on 9 and the coefficient 
of 9 must be separately equal to 0. This yields the following set of ordinary differential 
equations (ODEs): 

V^i' - V^i - H^At! + H 1 'A 1 = (12) 



AM - HM + - faVo' = Cr 



'01 



(n/V) 



A\ 



rA{) -A 1 



■Hi') 



(13) 
(14) 



— (ript ) - ipi 
r 



-Ho')' 
r 



— (rA 1 ')'-A 1 
r 



-Mo') 

r 



Ax'Vx - AxVx + ipxHx - Hxipx - W H x" - V H i = 



(15) 
(16) 



VxA ' + faH ' - AxV ' - H^' - r,rH^' - V H ' = (17) 
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if)[Ai - if)iA[ + t}(Ai + rAi") = (18) 

E z + ^Ax - ^ + rj(A ' + rA ")} = 0. (19) 

Here the prime stands for the derivative with respect to r. The above system consists of 
8 ODEs for 8 unknowns of one variable. Therefore the ansatz (|TTJ) is compatible with the 
original system of partial differential equations. In our stationary case V x E = 0, which 
together with Eq. ( |TT1) yields a uniform ^-component of the electric field. 

One can also notice that Eqs. fllil) , fll5l) , (|18D and ( |19|) are the same as for the purely 



two-dimensional (2D) case corresponding to setting H = V = 0. Thus, the 2D equations are 
a limiting case of our 2| D case and, what is even more important, they are decoupled from 
the rest equations of the system. Solutions of Eq. ([14]) , (|T5|) , (ITSf) and (|T9|) have already 



been presented in Ref. [12| and we can simply use these solutions for solving our more general 
problem. It is worth noticing that, for known A\, Aq', ip\ and ipo , the system of equations 
(0), (0), (|T6| ) and (|T7D is linear in the functions Hi, Hq, V% and Vq. 



IV. IDEAL SOLUTIONS 

In this section Eqs. (|12j)-(]l9D are analyzed in the limit of vanishing resistivity rj = 0. 
We start by considering the equipotential case, i.e. the case where the electric potential 
is constant along the z axis so that E z vanishes as well. In this case, according to (|10|). 
[if), A] =0 and flTID the relationship 

ip = a. A, a = constant (20) 

is valid. Eqs. (0), ([T5|), (0) and ([IJ) then yield 

A 1 = C7 1 lnr + C7 2 , ipi = aA u (21) 



A = C 3 r 2 + C 4 In r, % = aA , (22) 

where Ci, Ci, C3 and C4 are arbitrary constants. These solutions describe a magnetic 
configuration with a field-aligned flow in the (r, 9) plane. In this configuration there are 
one magnetic null point and one stagnation point, both are located at r* = [r c , — (2Csr c 2 + 
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Cy/Ci]. Here r c = exp(— C2/C1) is denoted as a critical radius. 

It is not difficult to see now that Eqs. ([12]) and (|T7| ) can be satisfied in the ideal limit by 

Vi = Cl A u H x = c 2 A u (23) 

where c\ and c 2 are other arbitrary constants. Therefore V\ and Hi also vanish together 
with A\ and ipi at r c . Then the evaluation of fll3| ) at r c requires C = 0. Substituting of fl23"D 
into ([1|) and (0) yields 

Vq = h^—?{A - «Vo) - -^-r(aA - ^ ) = ci(C 3 r 2 + C 4 lnr) + C 5 , 

or — 1 cr — 1 

#o = — ^(aA) - Vb) - ^t(^o - a^o) = c 2 (C 3 r 2 + C 4 lnr) + C 6 , (24) 
o; z — 1 a z — 1 



where C5 and Cg are arbitrary constants. 
Evaluating B z and f 2 at r* we obtain 



B z \ r=r , = -2c 2 C 3 r c 2 ( lnr c - i) - ^ (c 3 r c 2 + ^ ) + C 6 
^ 2 | r = r * = -2ciCW f lnr c - 1) - ^ (c 3 r c 2 + ^' j + C 5 



(25) 



which in general do not vanish. Therefore, contrary to the 2D case, in 2^D we generally 
have neither nulls nor stagnation points in the limit E z = rj = 0. 

The velocity and magnetic field lines are presented in Fig. [I]. They are superimposed 
on the distribution of the magnetic z component. The poloidal components (r, 9) and 
the nonuniform parts of the toroidal components (z) of our v and B are proportional to 
each other but generally with different coefficients of proportionality. This means that in 
our case the flow reshuffles the magnetic field lines inside each of the magnetic surfaces 
A(r, 9) = constant but it preserves the surfaces themselves. It can be noticed also that the 
distribution of magnetic flux can be imagined as generated by three sources with alternating 
polarities lying on the plane y = 0. 

Consider now how the above ideal equipotential solution is modified in the presence of a 
finite z component of the electric field. In two dimensions it has been shown in Ref. [| that 
a non-vanishing E z causes the appearance of a singularity at r = r c in the current density, 
vorticity and azimuthal components of the magnetic and velocity fields. Indeed, one class 
of solutions can be written as follows 

A l =C l \nr + C 2 , Vi = olA x , (26) 




Bz max 



Bzr 



FIG. 1: Magnetic (solid) and velocity (dashed) field lines superimposed on the distribution of the 
z component of the magnetic field shown in gray half-tones. The parameters used in the plot are 
E z = rj = 0, r c = 0.2, a = 2,d = 1.2, C 3 = 0.2, C 4 = -O.O8C3, C 2 = -Ciln0.2, C 6 = 1, c 2 = 0.8, 
d = 0.05. 



. , a E z r a b , , 

4/ = — + -r + — , Vo'= , . , 

or — 1 A\ a ar a z — 1 A\ 



1 E z r b 

h ar + -, 

r 



(27) 



where a and b are arbitrary constants. From here one can see that the magnetic flux piles 
up at the separatrix r = r c to produce there the above mentioned singularity. Coming back 
to our 2~ D problem, we notice that Eqs. (|12"D , (pT3|),(p!6D and ( |17D do not depend on E z . 
Therefore the expressions (p4|) are still applicable for E z 7^ 0, if one uses for t4 ' and ipo the 
expressions (^3). The solutions for V±, Hi, Vq and H$' are then given by 



V x = Cl A 



i- 



Hx = c 2 A 1: 



(28) 



Ho' 



ci E z r c 2 ( b 
H \ ar + - 



1 A x 



■st 1 ci / 6\ c 2 ,E 2 r 

r ° = «rnr^iv (29) 



These expressions show that, as in the 2D case, the presence of a non- vanishing E z leads to 
the appearance of singularities in the distributions of physical values. 

In particular, there is a singularity in the current density distribution (Fig. 0) as well as 
in the distributions of Aq', ipo', Hq' and Vq '. This means in turn that both the toroidal (z) 
and poloidal components of the current density become singular, which is also clear from 
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FIG. 2: Distribution of the current density J = |j| in the case of rj = and finite E z . The other used 
parameters in the plots are c\ = 0.5, C2 = 0.4, r c = 0.2, a = 0, b = 1, a = 8/9, C\ = —0.9/ In 0.2, 
C 2 = 1, d = 0.05. 

their explicit expressions: 

■ -Ei 

Jr j 

r 

jo = -H[6 - H' , (30) 
jz = — rA . 

In other words, these 2|D solutions inherit the singularities from the corresponding 2D 
solutions although the interpretation of the inherited singularities compared to the original 
ones differs in several important respects. 

As was shown above, the equipotential 2D magnetic configuration has a null point or, if 
one extends the system to three dimensions, a null line parallel to the invariant z direction. 
If an electric field E appears in the system, it may have just one z component due to the 
assumed two-dimensionality of the flow. The steadiness of the flow requires also that this 
electric field must be uniform. The latter implies in particular the presence of the electric 
field at the place where the null point (or null line) was initially located in the equipotential 
configuration. According to the frozen-in law condition, however, the finite value of E is 
sustained by the inductive field — v x B only. So, kinematically, the presence of E ^ at a 
magnetic null point would compel the velocity v to be infinite there. In our self-consistent 



S 



approach, incorporating both kinematics and dynamics of plasma, any velocity singularity 
may coexist only with an appropriate magnetic singularity, appearing at the same point 
in the result of the corresponding force balance. This is the reason why the magnetic null 
point of the equipotential configuration transforms into the proper magnetic singularity 
when passing in the considered family of exact 2D solutions to the configuration with a 
non-vanishing electric field. 

The situation becomes different for our 2^D configurations, which generally have no longer 
null points. Therefore we cannot appeal to the above "null-point" argument to explain 
the origin of the inherited singularities. The desirable explanation in fact can be found 
considering the consequences of the presence of field lines with a finite longitudinal voltage 
dropEl, denoted also as singular magnetic field /mes0. Our 2|D equipotential configuration 
has in place of the null line a straight magnetic field line parallel to z axis. Passing to 
non-equipotential configurations in our family of solutions leads to the appearance of a 
constant z component of the electric field, which inevitably creates the above mentioned 
voltage drop along a straight magnetic field line. Kinematically, such a voltage drop in 
a plasma with an infinite conductivity may be sustained by a suitable velocity singularity 
on l5 a Our self-content consider, incorporating plasma dynamics, shews that such 
a velocity singularity gives birth to an appropriate magnetic singularity by analogy with the 
2D case. Thus, our stationary 2|D ideal MHD solutions provide a particular but explicit 
realization of singular magnetic field lines, whose properties were kinematically described 
first by Schindler et al.0 in a more general non- stationary case, where the voltage drop was 
localized at a finite part of such lines. 

The above qualitative consideration shows that the appearance of singularities at the 
null line or singular magnetic field line can be anticipated if one combines the results of an 
analysis based on both the kinematics and the dynamics of a plasma. However, our explicit 
solutions reveal a much less obvious feature of this process, namely, that the singularity 
appears not only at the null or singular magnetic field lines but also at the whole magnetic 
separatrix surface r = r c containing such lines. In a more simple neutral X-type point 
configuration a similar fact follows from the frozen-in law and incompressibility conditions 
if one assumes also that the resulting plasma flow crosses one of the two separatrices and 
is parallel to the otherS. These conditions are sufficient for the appearance of a singularity 
along the separatrix which is not crossed by the flow. The self-consistent incorporation of 
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the plasma dynamics just makes the type of such a singularity more precise. One can also 
see this from the ideal MHD limit of Craig-Henton solution! and from the same limit of 
more general solutions!. It is not difficult to show that the considered point of view remains 
valid for our 2D and 2|D solutions in curvilinear geometry as well. 

Thus, the above discussion suggests that the "extension" of the singularity from the null or 
singular field line to the whole separatrix surface is somehow effected by the incompressibility 
condition. This suggestion, however, looks not convincing enough if one remembers that in 
the three-dimensional case the velocity singularity at the fan separatrix surfaces of magnetic 
nulls appears already in a purely kinematic approachS without any involvement of the 
incompressibility condition. So the complete clarification of the nature of the separatrix 
singularities is still a matter of the future development of the theory. Irrespective of the 
ultimate answer to this question, it is important to study how such singularities are resolved 
in the framework of a self-consistent MHD approach by a finite resistivity, which is an issue 
of the next section. 



V. RESISTIVE SOLUTIONS 



Let us now consider the complete system (|12[)— (|19[) . For the functions A%, ipi, A' Q and iJj' 
we can use the following solutions 



A x = Cilnr + C 2 , il>i = aA 1 , 



(31) 



A '(r) 



cxp 



-(a 2 -l)(Ci lnr+C 2 ) 2 
2rja 



TjT 



r E z t - Ai^- (at + b/t) dt 



exp 



-(a 2 -l)(Ci\nt+C 2 y 
2r/a 



(32) 



Mr) 



1 

a 



a n \ 1 - a f b 

A (r) (ar + - 

a \ r 



(33) 



derived in Ref. IT21. These solutions describe a 2D reconnective annihilation in a curved 



current layer formed on one of the magnetic separatrices when the other is crossed by a 
sheared flow. The magnetic and velocity fields have a null point and a stagnation point, 
respectively, whose positions, contrary to the case considered in Sec. [TV1 are in general not 
coincident. Considering now the 2^D problem, we first notice that Eqs. (p8|) are solutions 
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of the system for a finite resistivity as well. Therefore the problem is reduced to finding 
solutions for H' and Vq. By means of Eq. flT3| ) we can express V ' in the form 

Tl AxH'v-H^ + Vx^ 



V' 



■0i 



which, with the help of ( |31~1) and (|28|) , can be reduced to 



1 



V^ = -(H' -c 2 A' + cM. 



a 



(34) 



(35) 



Inserting this expression into Eq. (|I?D and again using (|3~ID and (p8|), we obtain the following 
equation 



-A 1 + r ] 



H'o + Ax (£i + C2 )^-(^ + c 1 )A( ) 



a 



c 2 



a 



0. 



(36) 



This linear equation can be solved by the method of the variation of parameters with the 
following result: 



ar/r 
exp 



— / da^iaAo'-M+H^Ao'-aik')) 
r)r J 



(« 2 ~ 1) 



Ax , . 
—at ) + c 5 



exp 



^1 f±i. 



(37) 



By means of (p8|) , (|3TD and (|3^) the solution (g^) can be written as follows 



Hq = exp 



(a 2 -!)^ In r + C 2 ) 
2?7Q; 



-c 2 - 



Q 



as 



(Cx In s + C 2 )exp 



ds I cx(a 2 — 1) I tpQ — as — 



+ c 5 



(a 2 -l)(Ciln5 + C72) 
2r/a 



(38) 



The resulting magnetic field configuration is represented in Fig. |^. The plot refers to 
the particular case where the line (r = r c ,8 = 0) is a magnetic null line. A more general 
configuration without this null line may be obtained by simply adding any constant value to 
the corresponding B z distribution. This would give us an example of the resistive solution 
which in the limit of the ideal MHD describes the above mentioned configuration with 
singular magnetic field line. 

In Fig. | the absolute value of the current density is plotted. Comparing this plot with 
the one shown in Fig. [|, one can see the effect of introducing a finite resistivity, which 
indeed resolves the singularity at the critical radius. 

11 




FIG. 3: Integral lines of the poloidal magnetic (in black) and velocity (in gray) fields, superim- 
posed on the distribution of the z component of the magnetic field shown in gray half-tones. The 
parameters used for the plot are rj = 10~ 2 , E z = 0.5, c\ = 0.5, C2 = 0.4, r c = 0.2, a = 0, b = 1, 
a = 8/9, Ci = -0.9/ In 0.2, C 2 = 1, d = 0.05. 




FIG. 4: Distribution of current density for the resistive case. The parameters used for this plot 
are the same used for the plot in fig. |3[ 

As already noticed in Ref. [| the configuration of the magnetic field in the xy plane described 
by our resistive solutions is of particular interest for modeling the reconnection process in a 
special large class of solar flares. One can see from Fig. ^| that the magnetic field here is gen- 
erated by three sources with alternating polarities on the plane y = 0. Observations^! show 
that a large fraction of solar flares occurs namely in configurations with three photospheric 
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magnetic sources. The 2^D solutions presented in this paper provide a generalization which 
makes the previous 2D model more realistic. 

VI. CONCLUSIONS 

We presented a class of steady solutions of resistive incompressible MHD equations in 
cylindrical coordinates. These solutions describe three-dimensional velocity and magnetic 
fields with a translational invariance along the z axis in a half-space corresponding to the 
solar corona. Both fields are represented as a superposition of poloidal and toroidal (z) 
components, so that the resulting fields have no null points in the volume. The poloidal 
fields, however, have one null line along which two separatrix surfaces intersect. One of 
the separatrices is the same for poloidal velocity and magnetic fields and it is a segment of 
cylinder with the edges at the photospheric boundary. The other separatrices are different 
and they intersect the cylindrical separatrix along the null lines of poloidal fields in such 
a way that these lines lie at some distance from each other. The resulting plasma flow 
intersects the non-cylindrical separatrix of the poloidal magnetic field to produce a strong 
shearing motion along the cylindrical separatrix, where a strong current layer is formed. 
Thus, our solution is a curvilinear analog of the planar reconnective magnetic annihilation 
considered earlier in Ref. § and generalized in Ref. 0. In addition to the curvature it 
has two more features of interest. First, there are no n- vanishing and non- uniform toroidal 
components of the velocity and magnetic field in the corresponding configuration. Secondly, 
the resulting distribution of magnetic fields on the photosphere represents three areas of 
alternating polarities. This feature is very interesting for modeling a wide class of solar 
flares which have three photospheric magnetic sources of different polarities. 
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